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ABSTRACT 

We present X-ray imaging, timing, and phase resolved spectroscopy of the 
anomalous X-ray pulsar IE 2259+586 using the Chandra X-ray Observatory. The 
spectrum is well described by a power law plus blackbody model with F = 3.6(1), 
kTBB= 0.412(6) keV, and Nn= 0.93(3) x 10^^ cm'^; we find no evidence for 
spectral features (0.5 — 7.0 keV). We derive a new, precise X-ray position for 
the source and determine its spin period, P = 6.978977(24) s. Time resolved 
X-ray spectra show no significant variation as a function of pulse phase. We 
have detected excess emission beyond 4" from the central source extending to 
beyond 100", due to the supernova remnant and possibly dust scattering from 
the interstellar medium. 

Subject headings: pulsarsdndividual (IE 2259+586) — stars:neutron — X-rays:stars 

1. Introduction 

X-ray pulsars, typically, have hard spectra and a broad range of observed spin peri- 
ods (2 ms - 3 hr). There is, however, a group of at least 5 X-ray pulsars (4U 0142+61, 
IE 1048.1-5937, IRXS J170849-40091, IE 1841-045, IE 2259+586) that have character- 
istics similar to one another, but distinct from those known to be accreting binary pulsars. 
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The existence of this class of pulsars, dubbed the 'Anomalous X-ray Pulsars' (AXPs), was 
first identified as such 6 years ago (Mereghetti & Stella 1995; van Paradijs et al. 1995). AXP 
characteristics include: (i) predominantly constant X-ray luminosities of order ~ 10^^ ergs 
s~^; (ii) relatively soft, two-component X-ray spectra (blackbody kTsB ~ 0.5 keV; power-law 
photon indices ~ 2.5—4); (iii) small Galactic scale height (-Zrms ~ 56 pc); (iv) association 
of 2 out of 5 (IE 1841—045, IE 2259+586) with young (~ 10"^~^ year) supernova remnants 
(SNRs) (Gaensler et al. 2001); (v) pulse periods within a narrow range (5 - 12 s); (vi) rapid, 
nearly constant spin-down rates (10^^^ — 10~^^ s s~^); and (vii) no clear evidence for a binary 
companion fiUing its Roche lobe (see e.g. Mereghetti 1999, for a review). 

IE 2259+586 was discovered by Gregory & Fahlman (1980) and lies along the line of 
sight of an X-ray and radio bright galactic supernova remnant G109. 1-1.0 (CTB 109) (Rho 
& Petre 1997; Hughes et al. 1984). There is no accurate distance measurement for the SNR; 
Sofue et al. (1983) and Hughes et al. (1984) have estimated distances of 4.1 and 5.6 kpc, 
respectively, using the surface brightness-diameter relation in different wave bands. The X- 
ray properties of IE 2259+586 (and the associated SNR G109. 1-1.0) have been investigated 
using BeppoSAX, ROSAT, BBXRT, and ASCA (Parmar et al. 1998; Rho & Pctrc 1997; 
Corbet et al. 1995). The pulsar spectrum has often been described as an absorbed powerlaw 
(PL) with a photon spectral index of F = 4.0 and a blackbody (BB) component with a 
characteristic temperature of A;Tbb= 0.43 keV with galactic absorption column density of 
A^H=0.9 X 10^2 cm-2 (Rho & Petre 1997). 

The best previous position for IE 2259+586 was a = 23^ 01"^ 08".44, 6 = 58° 52' 44". 1 
(J2000) with a 95% confidence error radius of 2". 2 established by HuUeman et al. (2000) 
using ROSAT HRI observations. They found no optical counterpart with the Keck, down 
to hmiting magnitudes oi R — 25.7 and / = 24.3. Coe et al. (1994) set an upper limit with 
the VLA to the 1.5 GHz flux of 50 //Jy. 

We present here an analysis of our Chandra observation of the source. In §2 we describe 
the observations and present imaging, timing and spectral results. The analysis and deriva- 
tion of a new precise position of a = 23'^ 01"^ 08^295, S = +58° 52' 44".45 (J2000) with a 
99% confidence error radius of 0".60, together with the K-band identification of IE 2259+586 
is discussed in the companion paper by HuUeman et al. (2001). 

2. Observations and Results 

IE 2259+586 was observed on 11 January 2000 using the Advanced CCD Imaging 
Spectrometer (ACTS). Data were collected sequentially in two different observing modes: 
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timed exposure (TE) mode (19 ks) and continuous clocking (CC) mode (12 ks). Data 
obtained in TE mode allow for two-dimensional imaging. Accurate spectroscopy of the 
bright pulsating target is limited due to pulse pile up. Pulse-phased spectroscopy is further 
limited by the 3.24 s time resolution. In CC mode the amount of pileup is negligible due to 
its 2.85 ms time resolution, allowing for both accurate time integrated and phase-resolved 
spectroscopy. Furthermore, one can exploit the one-dimensional image to search for extension 
of the central source. 

The source was positioned on the nominal target position of ACIS-S3, a back illuminated 
CCD on the spectroscopic array (ACIS-S) with good charge transfer efficiency and spectral 
resolution (Townsley et al. 2000). In addition to S3, four front illuminated CCDs were active 
(12, 13, S2 and S4, see also Figure 1). The focal plane operating temperature was —110° C. 

Standard processing of the data was performed by the Chandra X-ray Center. The 
data were filtered to exclude events with ASCA grades 1, 5, and 7, hot pixels, bad columns, 
and events on CCD node boundaries. We removed the serial clocking noise streaks from 
S4 using the destreak program. We examined the processed data and found no times with 
bad aspect. The S3 fight curve was inspected in a region offset from the AXP to identify 
periods of high background rates: we removed segments where the background exceeded a 
3(j threshold about the mean. The resulting useful observing times are 15600 s and 9000 s 
for TE and CC mode data, respectively. 



2.1. Extended Emission Search 

Searching for extended emission in the vicinity of IE 2259-1-586 was difficult. One must 
account for numerous contributions to the observed flux. These contributions include not 
only the bright source itself and the instrumental background, but also the SNR, stars in the 
field, and diffuse emission from the galactic plane. One must also account for the possibility 
of a dust scattering halo and/or an X-ray nebula. 

To avoid issues related to the effects of pileup in the TE mode data, we first utilized the 
one dimensional image from the CC mode. We generated a time-integrated image (0.5 — 7.0 
keV) minus any point sources other than IE 2259+586. We then subtracted an average coimt 
rate intended to remove instrumental and cosmic diffuse X-ray background, as measured by 
S3 (Markevitch & Vikhlinin 2001, and references therein). Next, we constructed what will 
be referred to as the "pulsed" image. We took the observed pulse profile, normalized it to 
a mean of zero and convolved it with the event fist of S3. Each count recorded on S3 was 
assigned a phase. Events with a phase near pulse maximum, regardless of their position 
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on the chip, received a positive weight and hkewise, those near pulse minimum, a negative 
weight. In doing so, we remove all emission components in our image that do not vary 
in phase (i.e., everything except the central pulsar). Note that the time delay for photons 
scattered by the interstellar dust is on average minutes, much longer than the pulsar period. 

To improve statistics, we folded the one-dimensional images about the common cen- 
troid and accumulated "quasi-radial" profiles. The timc-intcgrated (CCt) and pulsed (CCp) 
profiles are shown in Figure 2. We find that the pulsed radial profile is consistent with the 
MARX (Wise et al. 1997, v3.01) derived point spread function (PSF). 

Next, we collapsed the TE mode observation to mimic the 1-D CC image (Figure 2, 
TE). We find that the TE profile completely overlaps the CCt profile beyond ~ 4", marking 
the radius beyond which the effects of pile-up are negligible. In both profiles we find excess 
emission beyond ~ 4" which is likely due to the combination of the underlying SNR, a 
dust-scattering halo and potentially an X-ray plerion. We conclude that the majority of 
the IE 2259-1-586 X-ray flux is contained within a radius of ~ 4". Inspection of the spectra 
extracted from S3 in five annuli, centered on the pulsar position, and with inner and outer 
radii of 4" — 20", indicates the presence of continuum radiation in the vicinity of the pulsar. 
Detailed modeling of the angular and spectral distribution is underway for an accurate 
accounting of the light, to establish the relative contributions from the SNR, dust scattered 
light from the pulsar, and any possible plerionic emission. This work will be presented 
elsewhere. 



2.2. Timing Results 

High resolution timing with ACIS is only possible using CC-mode data. The CC mode 

event times denote when the event was read out of the frame store, not when it was detected. 
We corrected for this effect by assuming that all photons were originally detected at the 
nominal target position. We removed the variable time delay due to spacecraft dither and 
telescope flexure using the ACIS CC mode absolute time corrector^. The event arrival times 
were then corrected to the solar system barycenter using axbary which utilizes the JPL 
planetary ephemeris DE-200. 

The data were divided into ten ~ 1200 s intervals and individual pulse proflles were 
derived using epoch folding. These ten proflles were compared to the pulse proflle using all 
the data and the relative phases were flt with a linear function. The resulting pulse period 



^http: / /wwwastro. msfc.nasa.gov/xray/ACIS / cctime/ 
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of 6.978977(24) s is referred to epoch MJD 51555.0. This period is consistent with the spin 

history of the source (Gavriil & Kaspi 2001). The pulse profile is shown in Figure 3a; the 
0.5—7.0 keV background subtracted peak-to-peak pulse fraction, as defined in Ozel et al. 
(2001), is 35.8 ±1.4%. 

Previously, Mereghetti et al. (1998) have established limits as to the existence of a 
binary companion. They concluded that a putative companion must be a white dwarf or a 
He-burning star with M < O.SMq (for a companion star filling its Roche lobe). We have 
also searched for evidence of binarity. Binning the data in 50 s intervals, a total of 177 
pulse profiles were cross-correlated to the average pulse profile, and pulse phase offsets were 
determined. We then searched for the sinusoidal signature of a circular binary orbit, for 
periods ranging from 100 to 5 x 10^ s using the method described by Wilson et al. (1999). 
We found no evidence for binarity and set a 99% confidence limit of sin i < 70 It-ms for 
orbital periods in the range from 170 to 5000 s (see also Mereghetti et al. 1998). 



2.3. Spectral Analysis 

2.3.1. Phase Averaged Spectroscopy 

We use CC mode data to obtain an X-ray spectrum. We assume that the CC and TE 
spectral responses are identical. To test this assumption, we compared low count rate data 
from a particular region in both CC and TE mode and found reasonable agreement. We 
define the source by selecting an interval ±8 pixels (± ~ 4") (see also §2.1) around the 
peak fiux along the collapsed CC mode axis. The background was determined using two 
adjacent segments 12 pixels wide for a total area 1.5 times the source area; the background 
flux provides ~ 2% of the total flux in the source region. Source and background spectra and 
response files were generated using the CI AO (v2.1.1) tools dmextract, mkrmf, and mkarf. 
The extraction was performed in pulse invariant (PI) space (i.e., after the instrument gains 
were applied). The spectrum was grouped into bins that contained at least 25 events. 

AU spectral fits were limited to the 0.5 - 7.0 keV band and use XSPECvll.Ol (Ar- 
naud 1996), the photo-electric absorption coefficients of Balucihska-Church & McCammon 
(1992), and abundances of Anders & Grevesse (1989). We fit the data to a PL-I-BB model 

{x^/d.o.f. = 316.1/272); the fit is significantly improved {x^/d.o.f. = 262.5/258) by ig- 
noring the 1.9-2.1 keV band, which contains the iridium-edge structure in the telescope 
response. The best fit PL-I-BB parameters are F = 3.6 ±0.1, A';Tbb= 0.412 ± 0.006 keV, 
and Ny{_= 0.93 ± 0.03 x 10^^ cm~^, consistent with the same model results of Rho & Petre 
(1997). Assuming a source distance of 4 kpc, the 2-10 keV unabsorbed fiux and X-ray lumi- 
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nosity are 2.0 ± 0.2 x 10~^^ ergs cm~^ and 3.8 ± 0.4 x 10^^ ergs s~^ respectively. The 
PL component contributes 50% of the total 2-10 keV unabsorbed flux. Figure 4 shows the 
data, the best-flt model, and the residuals. We have also investigated single component PL, 
BB, and bremsstrahlung (BS) models and a combined PL+BS model and found that they 
are all statistically unacceptable at > 99.9, > 99.9, > 99.7, and > 96.4% confidence levels, 
respectively. 

We have searched for spectral lines. Features that appear below ~ 1 keV and in the Ir 
edge, are likely due to uncertainty in the spectral response. To set upper limits, we examined 
two small deviations at 0.7 and 5.0 kcV. The addition of a line, modeled by a Gaussian with 
an intrinsic width of 40 eV at 5.0 keV to the PL+BB model, results in an F statistic of 
3.3; the 90% confidence upper hmit on the fine fiux is 8.6 x 10~^^ ergs cm~^ s~^. Similarly, 
including a "fine" at 0.7 keV gives an F statistic of 1.3. We conclude that there are no 
significant line features. 

2.3.2. Phase Resolved Spectroscopy 

We fit the data from each of the ten pulse phase bins (Figure 3) with a PL+BB model 
parameters free to obtain their best fit values and with A'h constrained to be identical at 
each phase. The fitting gave a column density of 0.93 ± 0.04 x 10^^ cm~^ (and is identical 
to the time averaged value found in §2.3.1). For the remaining fits, A^h is constrained to 
be identical at each phase (i.e., linked). We next fit all phase bins with kT^B constrained 
to be identical at each phase and allow F and the normalizations to vary. The decrease in 
was not significant indicating the fits were statistically equivalent. However, there were 
significant variations in F (±0.2). The fits were repeated with F constrained with respect to 
phase and with /cTbb and normalizations allowed to vary. Again this resulted in a statistically 
equivalent fit but with significant variations of A;Tbb (±0.02 keV). The evolution of spectral 
parameters is shown as a function of pulse phase in Figures 3b and 3c. We fit the spectra 
with /cTbb and F constrained to be identical at each phase with only the normalizations 
free to vary with phase and find that the fits are not statistically equivalent. These results 
indicate that there is marginal (~ 3a) evidence for evolution in either F, /cTbb or some 
combination thereof with phase, but we are unable to determine which. Finally, we impose 
an additional constraint where the ratio of the PL and BB model normalizations (PLNorm, 
BBTVorm) are constant. Our fits show no statistical improvement in complicating the model 
by allowing the normalizations to be independent. Fit statistics for the models discussed are 
given in Table 1. 
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3. Discussion 

We have observed the anomalous X-ray pulsar, IE 2259+586 with Chandra and deter- 
mined a new, precise position of the source of a = 23'* 01™ 8''.295, 6 = +58° 52' 44". 45 
(J2000) with a 99% confidence error radius of 0".60. We report in the accompanying paper 
(HuUeman et al. 2001) the details of the derivation, which led to the identification of a near 
infrared (NIR) counterpart at the Chandra position. This is the second AXP counterpart 
to be found in NIR wavelengths (HuUeman et al. 2000). 

We find that of the spectral models we have considered, the PL+BB one is statistically 
preferred for the pulsar spectrum. In addition, the superb Chandra spatial resolution allows 
us to determine that there is emission extending up to 100" from the pulsar, but we are 
unable to determine its precise origin - whether from the supernova remnant, and/or a dust 
scattering halo, and/or a plerion. 

We acknowledge support from the following grants: MX-0101 (C.K.), NAG5-9350 (RW.), 
GO0-1018X (S.P.). We thank F. Paerels and J. Vink for many useful discussions. 
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Tablc 1. Phase Resolved Spectral Fit Results for IE 2259+586 



Model Settings' 
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*Model normalizations for each phase bin are free (unconstrained) to obtain its best 
fit value while iVn is linked (constrained to be identical at each phase) for all model 
variations. 

^The numbers in parenthesis correspond to the F— statistic needed to claim that the 
given model is significantly better at the 95% and 99% confidence levels, respectively. 
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Fig. 1. — Chandra ACIS 0.5 — 7.0 keV image processed as described in § 2.3. The bright 
hne that passes throught the source is the ACIS transfer (or trailed) image (Weisskopf et al. 
2000). The false grey scale represents the number of counts detected. The ROSAT image 
(Rho & Petre 1997) is superimposed for comparison and to give a better view of the SNR. 
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Fig. 2. — Radial surface brightness profiles of the total (CCt) and pulsed (CCp) emission 
(0.5—7.0 keV). The dashed lines are the simulated Chandra point spread function. Downward 
pointing arrows denote 2a upper limits to the count rate. The triangles are TE mode data. 
Note the agreement between the total CC profile and the TE profile at radii > 4". 
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Fig. 3. — Top panel (a): The pulse profile (0.5 — 7.0 keV). Middle panel (b): Variation of 
the blackbody temperature from fitting each phase bin with a PL+BB model. The time 
averaged value is denoted by the dashed line. Bottom panel(c) : Variation of the PL index. 
All errors shown denote la. 




Fig. 4. — The best fit spectral model (PL+BB) and residuals in units of a. The feature 
at ~ 2.0 keV appears to be due to a small shift of the location of the Ir absorption edge 
between the response and the data. 



